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Abstract 



The evaporation coefficients of water in air and nitrogen were 
found as a function of temperature, by studying the evaporation of 
QQ I pure water droplet. The droplet was levitated in an electrodynamic 

CN ■ trap placed in a climatic chamber maintaining atmospheric pressure. 

Droplet radius evolution and evaporation dynamics were studied with 
high precision by analyzing the angle-resolved light scattering Mie in- 



QQ ' terference patterns. A model of quasi-stationary droplet evolution, 

^D ■ accounting for the kinetic effects near the droplet surface was applied. 

In particular, the effect of thermal effusion (a short range analogue 
of thermal diffusion) was discussed and accounted for. The evapo- 
^ ■ ration coefficient a in air and in nitrogen were found equal, a was 

5t , found to decrease from ~ 0.18 to ~ 0.13 for the temperature range 

from 273.1 K to 293.1 K and follow the trend given by Arrhenius for- 
mula. The agreement with condensation coefficient values obtained 
with essentially different method by Li et al.^ was found excellent. 
The comparison of experimental conditions used in both methods re- 
vealed no dependence of evaporation/condensation coefficient upon 
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droplet charge nor ambient gas pressure within experimental param- 
eters range. The average value of thermal accommodation coefficient 
over the same temperature range was found to be 1 it 0.05. 

Keywords: Mie scattering, evaporation model, Arrhenius formula. 



1 Introduction 

Many problems of science and technology are related to the evaporation 
from droplets and condensation on them. Cloud and aerosol microphysics 
together with construction of climate models, ^"^ electrospraying,^ combus- 
tion,^ jet printing (compare^) and spray painting (compare^) are just some 
areas of relevance. Though they concern large sets of coexisting droplets, the 
understanding of transport processes at the surface of a single droplet is vital 
for solving them properly. Mass and heat transport processes at (nearly) flat 
interface can be efficiently modeled as a diffusion phenomenon. However, the 
evolution of small droplets is significantly infiuenced by effusion, which takes 
place in effectively collision-free region in the very vicinity of the interface 
(up to the mean free path of surrounding gas molecules). In order to ac- 
count for this phenomenon, a so called evaporation (condensation) or mass 
accommodation coefficient a is introduced besides the diffusion coefficient. 
Likewise the thermal conductivity coefficient should be accompanied by a 
thermal accommodation coefficient a^. These coefficients describe transport 
properties of the liquid-gas interface. The mass accommodation coefficient 
can be perceived as the probability that a molecule (e.g. water) impinging 
on the interface from the gaseous phase side enters into the liquid phase and 
does not rebound. Analogically, the thermal accommodation coefficient de- 
termines the probability that a molecule impinging on the interface attains 
thermal equilibrium with the medium on the opposite side. The considera- 
tions of evaporation and condensation coefficients are considered equivalent 
and the values of these coefficients - equal. ^ Both a and ay coefficients are 
phenomenological and should describe only the properties of the very inter- 
face. All other processes influencing mass and heat transport, such as chem- 
istry of the interface or the electrostatic interactions should be accounted 
for separately. ^° It is agreed, however, that a might possibly exhibit some 
temperature dependence. ^'^^ 

Many attempts have been made over nearly a century, to determine ex- 
perimentally the values of a and a-jn for water, but the results obtained by 



different authors spanned from ~ 0.001 to 1 for a and from ~ 0.5 to 1 for ut 
(see e.g.^'^^"^^ and^'^^'^° for revues). A variety of experimental methods was 
used. Both condensation on and evaporation from the surface of bulk liquid, 
liquid films, jets, and droplets were investigated in various environments (vac- 
uum, standard air, passive or reactive atmospheres) under various pressures 
and for various water vapor saturations. Small droplets, such as encountered 
in clouds, have been favored, since kinetic effects manifest strongly for them. 
Suspended droplets, trains of droplets, clouds of droplets and single trapped 
droplets were studied. 

We must admit that in our studies we have also experienced the flow of 
kinetic coefficients values in time. We have tried to overcome it. We will 
discuss the possible sources of the divergence of results in section 15. 1[ 

The measurement of temperature dependence of a was rarely attempted, 
since the large divergence of obtained a values obscures the effect. Two 
recent studies by Li et al.^ and by Winkler et al.^^ can serve as an ex- 
ample. The authors of the first study (Boston College /Aero dyne Research 
Inc. group) found that a decreases with temperature within the temperature 
range between 257 and 280 K. The authors of the second study (University of 
Vienna/University of Helsinki group) claim that a and ut exhibit no temper- 
ature dependence between 250 and 290 K. The comparison of these results 
can be found in.^° 

In this paper we present our new experimental results of evaporation 
coefficient of water in air, as well as our reprocessed results for water in ni- 
trogen (compare^^), versus temperature; both under atmospheric pressure. 
The results for air and for nitrogen are consistent, which also indicates that 
the presence of small amounts of such soluble and/or reactant gases as CO2 
in the ambient air, does not influence the value of kinetic coefficients. In 
comparison to our previous works we refined our data processing which en- 
abled us to determine the droplet radius with higher accuracy and trace its 
evolution with higher confidence. Smoother radius derivatives enabled us, in 
turn, to employ direct fitting of the model in finding the kinetic coefficients 
and avoid most approximations. We also operated on a larger set of exper- 
imental runs. This yielded correction of the values of kinetic coefficients we 
obtained previously, and revealed different temperature dependence. These 
results turned out to be in excellent agreement with the results of BC/ARI 
group: the values of a coincide within the temperature range from ~ 273 
to ~ 281 K and our results extending towards higher temperature follow 
the temperature dependence found by BC/ARI group. Since BC/ARI group 



and our results together span over larger temperature range, the accuracy of 
finding temperature dependence of a could be improved. 

2 Experiment 



PC with multi-channel A/D & 

D/Acard, IEEE-488 card 

and frame-grabber 



piezo injector driver 



>U^ 



droplet injector 
gas outlet from top 



632 nm laser beam 



shutter 



Peltier elements 
driver 



_ temperature and 
humidity sensors 




polarizers 



r^ 



sine-wave generator 
and HV amplifier 



DC power supply 



electrodynamic trap 
lens 

532 nm laser beam 



shutter 



gas inlet from below 

water-cooled Peltier elements 
objective lenses 



video recorder 



gain control 



Figure 1: Experimental setup - top view. 

The experimental setup is presented in figure [H It consisted of a hyper- 
boloidal electrodynamic quadrupole trap (see e.g.^^), kept in a small climatic 
chamber. The high resistivity electric circuit of the trap drive enabled us to 
operate in a humid atmosphere. A detailed description of this apparatus can 
be found in^^'^^ and of further modifications in.^^'^^ 

The droplets were introduced into the trap with a piezo injector, similar 
to constructed e.g. by Lee et al.^^ or Zoltan.^^ The injection timing was 



controlled relative to the trap driving AC signal. By choosing the proper 
injection phase, the sign and, to a certain extent, also the value of the charge 
of the injected droplet could be controlled. The initial temperature of the 
droplet was that of the chamber. 

In each experiment, the chamber was first flushed with clean, dry nitro- 
gen, and then filled with a mixture of nitrogen/air and water vapor. The 
temperature in the chamber was monitored and stabilized. A zone type 
temperature control enabled us to eliminate vertical temperature gradients. 
Horizontal gradients were found negligible. 

The humidity in the chamber, but outside of the trap was monitored con- 
tinuously with semiconductor sensors. Due to poor vapor exchange through 
trap openings accompanied by injecting liquid water into the trap, the humid- 
ity inside the trap could not be inferred directly from those measurements. 
The value of humidity in the trap found as fitting parameter (see section H]) 
turned out to be higher by several percent than sensors readings. Resorting 
to fitting method was inevitable, since the humidity accuracy required for 
the correct assessment of kinetic coefficients was inaccessible via any on-line 
sensor measurement. On the contrary, analyzing droplet radius evolution 
seems a highly accurate method of measuring relative humidity, surpassing 
any on-line methods. 

In our experiments we used ultra-pure water. The details about its initial 
parameters and sample preparation can be found in,^^ where we discussed 
also the absorption of impurities by ultra-pure water, and their influence 
upon the experimental results there. 

Droplet evolutions were studied with time resolved static light scattering, 
with green or red laser light. We found no inconsistency between the results 
obtained for both and we infer that the light wavelength had no influence 
upon the results. 

3 Evaporation model 

In order to interpret the experimental results, a model of evaporation was 
necessary. The model of evaporation we used was based on a generally ac- 
cepted model which can be found in such textbooks as.^'^*^'^^ It was a slightly 
rephrased and numerically reexamined version of what we had used previ- 
ously.^^''^° Below we discuss the details of the model equations we used, since 
the results may depend significantly on the apparently minute approxima- 



tions made. We also point to a certain approximation typically made that 
we found weighing heavily upon the results. 

The quasi-stationary evaporation of a free, motionless droplet larger than 
the mean free path of vapor molecules, can be easily described with the 
diffusion equations with boundary conditions defined by the thermodynamic 
conditions in the reservoir (far from the droplet). This part of the model does 
not rise much difficulties as long as the characteristic times of the process 
justify the quasi-stationary approach. ^^'^^ 

For small droplets, of the size comparable with the mean free path of 
vapor molecules, the language of diffusion is not appropriate. The transport 
of mass and heat below the mean free path distance from the surface must 
be perceived as effusion or evaporation into vacuum. The net effusive flow of 
vapor can be expressed as difference between outgoing and incoming effusive 
flows:^'^^ 

J = na^a [pe{r = a)v{Ta) - p{r = a + A)v(Ta+A)] , (1) 



where v{T) = JSRT/ttM is the mean absolute thermal velocity of vapor 
molecules for the temperature T; Ta is the temperature of the droplet (sur- 
face), Ta+A is the temperature of vapor at the distance A (comparable with 
the mean free path of the vapor molecule) from the surface. p{r) is the vapor 
density at the distance r from the droplet center while pe{r = a) is the vapor 
density at the droplet surface for the equilibrium conditions (steady state, 
no net flow). 

The usual approximation made is Ta = Ta+A (see e.g.^). It implies ne- 
glecting the slowing down of the mass transport by thermal effusion (a short 
range analogue of thermal diffusion). It should also be noted, that lifting 
the temperature dependence of v introduces some additional temperature 
dependence into a. Unfortunately, discarding this usual approximation ex- 
cludes using standard solutions and substantially complicates calculations. 
To overcome such difficulties, we decided to introduce a simple correction 
of a at the end. We shall address this issue in detail later. Following the 
standard route, we compare effusive and diffusive flows aX r = a + A. Since 
these flows are equal and both are proportional to the vapor density gradient 
it is possible to write a compact expression: 
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is the Kelvin equation, accounting for the modification of equihbrium vapor 
density near the droplet surface due to the surface curvature and charge 
effects, ^^ and 
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is the Clausius-Clapeyron equation. The effective diffusion coefficient Dk 
accounts for the effect of effusion: 
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D is the diffusion constant for water vapor in nitrogen/air, Tr is the temper- 
ature of the reservoir, Q is the droplet charge, p^o and Pa are the equilibrium 
(saturated) vapor pressure above the flat interface and above the interface 
of the curvature radius a at a given temperature. S is relative humidity. 7, 
PLi M and q are the surface tension, density, molecular mass and the latent 
heat of evaporation of liquid water, e^ is the permittivity of vacuum, R is 
the universal gas constant. A^ deflnes the effective range of the gas kinetic 
effects. It is comparable with the mean free path of particles of surrounding 
gaseous medium A^. We assumed Ac = 4Aa/3.^ 

The change of droplet mass by evaporation (condensation) is associated 
with heat absorption (release), which manifests as temperature drop (rise) 
toward the droplet. The equation for the transport of heat can be presented 
in a convenient form: 
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(8) 



is the effective thermal conductivity of moist nitrogen (air) and A, ptv, ^n 
and cp are thermal conductivity, density, molecular mass and specific heat 
capacity under constant pressure of moist air /nitrogen respectively. Ay plays 
role analogous to Ac and was assumed as At = Ac + 4A/ (wcppTv)- Since 
in the vicinity of standard temperature and pressure, the partial pressure 
of water vapor can be neglected in comparison to that of air /nitrogen, it 
can be assumed that the heat is conducted to the droplet mostly by the 
molecules of air /nitrogen. In consequence, the flux of mass can be considered 
independently of the flux of heat and Ac associated with the transport of 
mass should be distinguished from At associated with the transport of heat. 
The direct influence of the droplet charge, through charge-dipole inter- 
action, upon the mass (or heat) transport was estimated to be negligible 
for droplet charge ensuring Coulomb stability (compare^''). Similarly, field 
emission did not take place for surface charge densities encountered in our 
experiments (see e.g.^^). The Coulomb explosion of the droplet is a threshold 
process and does not need accounting in the transport equations. 

4 Experimental data processing 

The procedure of the numerical processing of experimental data, which we 
found the most stable and yielding the most consistent results, relies on 
the direct fit of the model equations to the experimentally obtained droplet 
radius change rate a = da/dt as a function of droplet radius a. The data 
preparation procedure is presented below. 

The running radius of the droplet ai(tj) was obtained (off line) from the 
angularly resolved Mie scattering pattern for the time t, with the help of 
a gradientless fitting procedure ("library method"). Each droplet evolution 
yielded a sequence of a few hundreds data points indexed with i (see figure [2]). 
We had found that in order to obtain reliable results, significant care must 
be taken to ensure a high signal to noise ratio of the measurement. There 
happen data points misplaced to incorrect "evolution branch", associated 
with the Mie resonances that could not be handled with the method used 
(see description of the method^^). The accuracy of a single value of droplet 
radius Oj (except for misplaced points) was estimated as ±15 nm. The aj(ti) 
sequence was stripped to main "evolution branch" (indicated with arrow in 
figure [2]) and interpolated in order to obtain regularly spaced data points. 
The time derivative di(ti) was calculated (figure [2]). The aj(tj) evolution was 
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Figure 2: An example of temporal droplet radius evolution, before and after 
processing - top and bottom curve respectively. Derivative calculated from 
processed data - middle curve. A^ = 395, Tr = 283 K, patm = 1006 liPa, 



smoothed with low pass FFT filter and combined with the derivative in order 
to obtain di{ai). Finally di{ai) was smoothed (figure [3]). 

On the other hand, subtracting equation [2] from equation [7] leads to an 
equation binding Ta and a. For every experimental Oj this equation can be 
unambiguously numerically solved for Ta, yielding Ta(aj). This, on insertion 
into equation [2] yields at every experimental data point a numerically solv- 
able equation binding a and a^. Thus, a model prediction of a{ai) could be 
obtained. 

In order to find a, ot, S and Q, we minimized the function 
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using a gradient method. A^ was the total number of experimental data points 
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Figure 3: Droplet radius temporal derivative versus droplet radius, corre- 
sponding to figure [2|, before and after filtering (points). The result of model 
fitting - solid line. Fitting parameters: Sfu = 0.9762, Q = 3.7 x 10® elemen- 
tary charge units, a = 0.155, ax = 1- 



of the evolution, and xo was an arbitrarily chosen normalizing factor, a and 
S were found to be the essential parameters and could be unambiguously 
determined, while a^ and Q could be determined only with limited confi- 
dence. Since a and S had seemed partially interconnected, the minimization 
was performed very carefully, starting from various combinations of a and S 
{a larger, S smaller versus a smaller S larger) and accepted only if leading 
to the same results. The less relevant parameters were initialized as follows: 
aT = 1, (values above 1 were allowed; compare^^), and Q = SnJeo'yaf, where 
Qi corresponded to the smallest droplet radius observed in the evolution (no 
Coulomb instabilities during evolution). The resulting Q was much approx- 
imate, and we couldn't detect the eventual droplet charge loss (see eg.^^) 
by analyzing the evolution of the droplet radius. The minimization was also 
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hardly sensitive to ay, however a value close to unity could be inferred. Since 
for larger droplets (a > 6 /im) the kinetic effects as well as the effect of the 
droplet charge were negligible, only S was fitted in this range, as a first step, 
and then the minimization was extended towards smaller radii with a added 
as a parameter. Finally ot and Q parameters were added. The whole proce- 
dure exhibited best stability for 5* > 95%, since the evaporation was slower 
then (compare equation [2]) and thus: (i) the evolution of the droplet radius 
could be determined with high precision and (ii) the temperature jump at 
the interface AT was so small (compare equation [7j) that the model equations 
used were exact enough. It would be valuable to validate the procedure of 
finding kinetic coefficients using other liquids (such as ethylene glycol). Un- 
fortunately, the parameters such as diffusion constant are usually not known 
with adequate precision. On the other hand, after slight modification of the 
procedure, it should be possible to look just for diffusion constant, which we 
intend to do soon. 

4.1 Correction of a 

In order to estimate the influence of AT upon the obtained value of a, we 
apply an approximation T^+a = Tr to formula [H which is opposite to usu- 
ally apphed AT = 0, and we compare the results of both approximations. 
The approximation that we introduce means that we account only for ther- 
mal effusion while neglecting thermal diffusion. Since for our experimental 
conditions the temperature gradient was highest in the very vicinity of the 
interface (see^^), our approximation was legitimate. For simplicity we also 
assumed that the shape of distribution of vapor density was spatially con- 
stant and temperature independent. It implied p(r = a + A) = Sa+APiTn), 
where Sa+A = const represented relative humidity at r = a + A. If we require 
that the effusive flows calculated with each of the approximations are equal, 
we have: 
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Introducing Ta{ai) (see section H]) into equation [101 we can flnd a correction 
of a, where Sa+A is a (scaling) parameter. It is initiated as Sa+A = S and 
optimized so that a/a (AT = 0) — *> 1 for AT — > (larger T^ in case of our 
experiment; see the inset in flgureH]). The results presented in flgure [5] are 
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already corrected. In our case a significant (by a factor of nearly 2) correction 
was near the freezing point and by several percent at 276.5 K. The equation 
[To] is essentially approximate and leads to underestimation of a. It can be 
seen in figure [5] - our data points seem to lie slightly below the trend line. 
It turns out that for many reasonable experimental conditions the correction 
factor can be higher than 2. We shall discuss a few examples in section [3?T1 
Considering the approximations made, we estimate that for thermodynamic 
conditions encountered in atmosphere the accuracy of the correction factor 
should not be worse than several percent. 

5 Results and Discussion 

The raw results are presented in figure [Has a function of the droplet (surface) 
temperature. The kinetic coefficients should be presented as a function of the 
droplet (surface) temperature, since in general, due to evaporative cooling, 
it may differ significantly from the temperature of the reservoir. In case of 
BC/ARI group experiments, Tr — T^ < 2 K.^^ Though, in our case T^ — Tq, < 
0.7 K only, it is sufficient that some of our results correspond to supercooled 
water as well. 

The kinetic coefficients found for water droplets in nitrogen and in air 
were mutually compatible (see figure S]). It implies, that the gases absorbed 
by water from the air had negligible impact upon our measurements and 
generally there is no strong dependence upon the composition of the ambient 
atmosphere. 

The final results are shown in figure [51 There are values of evapora- 
tion coefficient we obtained (solid circles) and values obtained by BC/ARI 
group, taken from^ (hollow circles). The values of thermal accommodation 
coefficient we obtained are also presented (solid squares). Data points corre- 
sponding to our results were obtained by averaging the raw results (compare 
figure [1]). We also followed BC/ARI group and used the formula they de- 
rived basing on Transition State Theory (TST) (e.g. equation 7 in^^). Such 
formulation enables expressing the results in the language of thermodynamic 

potentials: 

ex 

exp{AGobs) , (11) 



1 — a 

where AGobs is the Gibbs free energy and its temperature dependence can be 
expressed as AGo6s = AHobs—TASobs- AHobs and TASgbs are treated just as 
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Figure 4: Non-averaged experimentally obtained values of a as a function of 
droplet surface temperature. Solid and open circles represent results obtained 
for nitrogen and air respectively. The corresponding calculated evaporation 
coefficient correction factors, due to the thermal effusion, are presented in 
the inset. 

parameters, their physical meaning is not clear (see discussion below). This 
formula is derived on an assumption, well justified with elegant experiments 
by Nathanson et al. described in,^* that the particles from the gaseous 
phase enter liquid via an intermediate surface state. Dashed line in figure 
|5] represents the fit we made to the results of BC/ARI group and our data 
points together. It yielded AHots = 4830 ± 150 cal/mol and ASobs = 20.3 ± 
0.5 cal/mol, which is within the limits of error equal to the values given in,^ 
i.e. AHobs = 4.8±0.5 kcal/mol and ASobs = 20.3±1.8 cal/mol. The accuracy 
of our fit (and so of the values obtained) is higher due to the larger number 
of data points. 

The comparison of our results with those of BC/ARI group indicates 
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Figure 5: Collected a and ax values as a function of droplet surface tem- 
perature. Solid circles and triangles: corrected evaporation coefficient and 
thermal accommodation coefficient respectively, obtained from our measure- 
ments; hollow circles: condensation coefficient measured by BC/ARI group. ^ 
Dashed line represents the ffi of the equation [TT] to the results of BC/ARI 
group and our data together. 

also that there was no perceivable influence of droplet charge upon kinetic 
coefficients. Vibrating orifice injector generates, at least in average, neu- 
tral droplets, while in our experiments with evaporating charged droplets 
it could be assumed that the charge was approaching its maximum value - 
the Rayleigh limit. Similarly, the comparison of aforementioned experiments 
reveals no measurable influence of ambient atmosphere pressure upon the 
value of kinetic coefficients. 

The temperature dependence of a, though obtained with essentially dif- 
ferent method, coincides with the results of BC/ARI group (see eg.^°). Our 
result extends into higher temperature range. Furthermore, we measured 
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evaporation coefficient while BC/ARI group measured condensation coeffi- 
cient. It supports the notion of equivalence of these coefficients. 

The thermal accommodation coefficient we obtained ut = 1±0.05 (figure 
|5]) agrees with both BC/ARI and UV/UH groups' results. However it is hard 
to asses the real uncertainty of ax] the statistical error we found may be too 
small (see section H]). Thus it is not possible to derive information on its 
temperature dependence. Recently, there seems to arise a general consensus 
that ut is close to 1, which means, that all the particles striking the interface 
thermalize. 

5.1 An attempt of results coordination 

Since it is quite improbable that all the kinetic coefficients measurements 
performed over the years were loaded with random errors, it must be assumed 
that the experiments, though accurate by themselves, measured different 
quantities. Many authors have tried to coordinate the results by pointing 
out what was really measured (see e.g.:^'^^). However there is no consensus. 
We shall also try to address this issue. 

The divergence of results obtained by different authors has been usually 
attributed to: (i) difficulties in accounting for various physical and chemical 
interfacial processes; (ii) effects of impurities, and especially surface active 
agents;^^ (iii) structure of the interface (dynamic surface tension, reaching 
the balance by the interface) and (iv) dependence of the coefficient value 
upon the model used (indirectness of measurement). It has been pointed 
out^'^^ that two classes of experiments could be distinguished: (i) with a 
quasi-static interface, yielding a < 0.1 and (ii) with a continuously renewing 
surface, yielding a > 0.1. However, such categorization requires defining the 
time scale. Such scale has not been agreed yet, neither the leading mecha- 
nism responsible for interface aging. For example, the characteristic times 
used in Molecular Dynamics (MD) studies are only hundreds of ps. This 
falls into a non-stationary interval, when the transients in the temperature 
and vapor density fields are starting to form. The Transition State Theory 
(TST) considerations of Nagayama et al.,^° seem to be in agreement with 
MD calculations and predict a ~ 1 around room temperature. However, it is 
worth noting, that, for example, stationary values of the surface tension are 
reached within milliseconds^^ and all these time scales are far below the char- 
acteristic timescale of cloud droplet growth process, which lie in the range of 
seconds (or even minutes).'*^ 
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Recently, Fukuta and Myers^'^ have noticed, that accounting for the ef- 
fect of moving gas-hquid interface ("moving boundary effect") can change 
the resuhing value of kinetic coefficients by several percent. In their work 
they managed to account for this effect in an elegant way. Though the ther- 
modynamical conditions and the velocity of the interface in our experiment 
were similar to theirs, in present work, we have decided to neglect the moving 
boundary effect, since the correction of mass accommodation coefficient we 
propose is much larger. 

In this paper we would like to point to a mechanism which falls within 
the 4th category - model dependent, however it is related to the issue of the 
characteristic timescale of the process and its distance from the thermody- 
namic equilibrium. Usually, authors are careful to estimate the characteristic 
times of mass and heat transport processes involved, in order to assure the 
proper description. It seems, that in some cases this alone can be somewhat 
misleading, because of the thermal effusion which we already mentioned. We 
shall consider four examples. 

In case of BC/ARI group experiments, the vapor-liquid contact lasts sev- 
eral milliseconds but the droplet is essentially in equilibrium with the reser- 
voir. In order to achieve temperatures below 273 K the evaporative cooling 
was used which inevitably caused temperature jump near the surface (up to 
2 K) and thermal effusion as a consequence. However, since the value of a 
was not obtained from the evolution of droplet radius, its value should be 
safe and no correction is needed. 

In our case, we selected for the analysis the droplet evolutions which 
lasted a few seconds which guaranteed that the process had been quasi- 
stationary in the diffusion time scale. For faster evolutions the temperature 
jump approached 1 K, and since a was obtained from the evolution, it had 
to be corrected by means presented above. 

In case of the experiment of UV/UH group, ^^ the evolution lasted ~ 50 
ms, which is shorter than in our case, but for the thermodynamic conditions 
they had, the process still could be regarded as quasi- stationary. However, 
the temperature jump of ~ 3 K could be expected for such evolution. This 
alone would require a correction of a by a factor of 2. Further overestimation 
might be caused by uncertainty of water vapor saturation. There are also 
rather few data points lying on a relatively flat curve, which as we know from 
our experience, causes the increase of measurement uncertainty. 

Lastly, in case of very interesting Fukuta and Myers experiment^'^ the 
evolution (condensation) lasted ~ 3 s (similarly as in our experiment). Since 
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the final droplet radius was ~ 2 /im, it can be inferred that a ^ 1 /im/s, 
which in turn yields temperature jump of only ~ 0.2 K. However, since the 
mass transport was relatively slow (supersaturation used was very small), 
the effect of even small temperature jump at the interface could be relatively 
large. According to our estimation (see expression [TOl) the correction of mass 
accommodation coefficient should be as high as 5! This would bring Fukuta 
and Myers result for NaCl and (NH4)2S04 at 277 K to a ~ 0.2, which agrees 
within the limits of error with ours and BC/ARI group results, even allowing 
for moving boundary effect which we neglected. 

6 Conclusions 

We conclude that it is feasible to obtain reliable values of evaporation coef- 
ficient by analyzing the evaporation of a small droplet. It requires however 
several tens of data points per evolution and droplet radius measurement ac- 
curacy of several nanometers. Generally accepted model of quasi-stationary 
evaporation seems sufficient for experimental data analysis in most cases. We 
found however that when evaporative cooling of the droplet becomes of the 
order of 1 K, it is necessary to consider the effect of thermal effusion, which is 
a short distance analogue of thermal diffusion. The kinetic coefficients found 
for water droplets in nitrogen and in air were mutually compatible. The evap- 
oration coefficient for the temperature range from 293.1 K down to 273.1 K 
was found to increase from ~ 0.13 to ~ 0.18 and follow the trend given by 
Arrhenius formula (see fTT!) with the parameters AHobs = 4830 ± 150 cal/mol 
and ASobs = 20.3 ± 0.5 cal/mol. This temperature dependence is in excel- 
lent agreement with the results of BC/ARI group, which concern condensa- 
tion coefficient, were obtained with essentially different technique for much 
lower ambient gas pressure and extend toward lower temperatures. The com- 
parison with BC/ARI group experiments enables to draw a few additional 
conclusions: (i) the evaporation and condensation coefficients are essentially 
equivalent; (ii) there was no measurable influence of ambient atmosphere 
pressure upon the value of kinetic coefficients in the range from ~ kPa to 
~ 100 kPa; (iii) there was no measurable influence of droplet charge upon 
the value of kinetic coefficients up to the Rayleigh stability limit. The value 
of thermal accommodation coefficient we obtained q;t = 1 ± 0.05 agrees well 
with recent results of many authors. 
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